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I. INTRODUCTION 

Since the birth of stereochemistry 150 years ago with 
Pasteur's discovery of handedness in molecules (Pas- 
teur, 1848) interest in chiral molecules has continued 
unabated. The term chirality was first coined by Lord 
Kelvin (Thomson, 1893) : 

"I call any geometrical figure, or group of 
points, chiral, and say it has chirality if its 
image in a plane mirror, ideally realized, can- 
not be brought to coincide with itself." 

Chirality permeates the entire fabric of the biological 
world. Indeed, life as we know it could not exist without 
chirality. The function of fundamental components of the 
cell, like actin, myosin, proteins, and lipids, relies upon 
their being chiral. The handedness of a molecule can af- 
fect its odor, potency, and toxicity. Thus the synthesis 
of a single enantiomer of a compound is crucial for the 
delivery of safe and effective pharmaceuticals and food 
additives. 



Since chirality is the absence of inversion symmetry, a 
structure is either chiral or it is not. However, just as 
the degree of order of a ferromagnet, which is either or- 
dered or not ordered, can be quantified, so the chirality 
of a structure can also be. A major theme of this article 
is the development of quantitative measures of chirality 
and its impact on physically measurable properties of 
materials. There is an agreed upon convention to answer 
the "yes" -"no" question of whether a molecule is chiral 
or not by identifying chiral carbons (to be discussed in 
more detail below) to which a handedness can be assigned 
via the Cahn-Ingold-Prelog rule (McMurry, 1992), which 
orders the chemical groups attached to the carbon ac- 
cording to their molecular weight. The identification of 
the handedness in this way, however, gives no indication 
of the magnitude, or even the sign of the optical rota- 
tory power this molecule will exhibit in solution. Indeed 
the magnitude and sign of the rotatory power depends 
on wavelength (though it is usually quoted in handbooks 
for the D-line of Na). Rotatory power provides only one 
quantitative measure of chirality, which could be used, 
for instance to assess whether one substance is more or 
less chiral than another. Similarly when a liquid crystal 
forms a cholesteric it is obviously chiral. It can be more 
or less chiral depending on whether the pitch is shorter 
or longer, respectively. There are many other such quan- 
titative measures. In this paper, we will define quantita- 
tive indices of chirality based on molecular geometry and 
show how they enter the determination of a particular 
observable, the pitch of a cholesteric liquid crystal. 

Microscopic chiral constituents have a profound effect 
on the macroscopic structures they form, striking exam- 
ples of which are common in liquid crystals (de Gennes 
and Prost, 1993). The simplest liquid crystalline phase 
is the nematic phase, characterized by long-range uni- 
axial orientational order of anisotropic molecules called 
nematogens, as shown in Fig. la. The centers of mass of 
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the constituent molecules are distributed homogeneously 
as in an isotropic fluid, but one of their anisotropy axes 
aligns, on average, along a common unit vector n called 
the director. Strongly biaxial molecules (Fig. 2) can in 
principle condense into a biaxial rather than a uniax- 
ial nematic phase with long-range biaxial oricntational 
order (de Gennes and Prost, 1993). In this phase, one 
molecular axis aligns along n, and a second orthogonal 
axis aligns on average along a second vector e perpen- 
dicular to n as shown schematically in Fig. lb. Biax- 
ial molecules can also condense into a uniaxial nematic 
phase with short-ranged biaxial correlations rather than 
long-range biaxial order as depicted in Fig. lc. 




FIG. 1. Schematic representation of (a) a nematic liquid crys- 
tal in which long molecular axes align on average along a spa- 
tially uniform director n. The nematogens of this phase can 
either be uniaxial (like Fig. 2a) or biaxial (like 2b). In the lat- 
ter case, the long axis is the c axis of length 2c. (b) Schematic 
representation of the plane perpendicular to n in a biaxial 
nematic. The b axes of nematogens align on average along e 
perpendicular to n. (c) Schematic representation of the plane 
perpendicular to n in a uniaxial nematic composed of biax- 
ial molecules. There is no long-range biaxial order but there 
are short-range orientational correlations that persist out to 
a correlation length £. The angle <f> measures the orientation 
of the the "biaxial"-axis of each molecule with respect to the 
a;- axis. 

If these nematogens are chiral or if chiral molecules 
are added to an achiral uniaxial nematic, the director 
n will twist creating the simplest twisted phase: the 
cholesteric or twisted nematic phase, the first liquid- 
crystalline phase to be discovered (Reinitzer, 1888). This 
phase is depicted schematically in Fig. 3. The director 
at position x = (ar, y, z) rotates in a helical fashion: 

n(x) = {cos qz, singz, 0) . (1) 

In equilibrium the twist wavenumber q assumes a pre- 
ferred value qo, which corresponds to a pitch P = 2Tr/q . 
Typically, the pitch can vary from hundreds of nanome- 
ters to many microns or more, depending on the system. 
Cholesterics with pitches on the order of 500 nm Bragg 
scatter visible light and appear iridescent. If chirality is 
added to a biaxial nematic, a cholesteric structure similar 
to that shown in Fig. 3 results, though the local molecu- 



lar order will be strongly biaxial. Other liquid-crystalline 
phases with macroscopic chiral structure include the blue 
phase, in which the director twists in all directions to 
produce a three-dimensional periodic crystal, and the 
smectic-C* phase. The dipole moments of the molecules 
in the latter phase become ordered as a consequence of 
their chirality and make this phase ferroelectric. Tech- 
nologies based on these ferroelectric liquid crystals show 
great promise for fast-switching, high-resolution displays. 




(a) (b) (c) 

FIG. 2. Representations of (a) Linear, (b) biaxial planar, and 
(c) chiral molecules. The molecule in (c) is a twisted "H" 
obtained from (b) by twisting about the long molecular axis 
(c-axis). It is both nonuniaxial and nonplanar as required for 
a chiral molecule. 
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FIG. 3. Schematic representation of a cholesteric liquid crys- 
tal, showing the helical twisting pattern of the local director 
n along the pitch axis (z in this case). The director rotates 
by 7r in half a pitch P/2. 

Models of chiral molecules are shown in Fig. 4. Of 
particular importance in chemistry is the tetrahedrally 
coordinated molecule, shown schematically in Fig. 4a, 
consisting of a central carbon atom with each of its four 
bonds connected to a different chemical unit. If any two 
of these chemical units are equivalent then the molecule 
has a mirror plane and is not chiral. If all the chemical 
units are different then the molecule is chiral and the cen- 
tral carbon atom is referred to as a chiral center. More 
complex molecules may have many such chiral centers. If 
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all atoms of a molecule lie in a single plane, that plane is a 
mirror plane, and the molecule is achiral; therefore chiral 
molecules must be three-dimensional. The converse does 
not hold: not all three-dimensional molecules are chi- 
ral. For instance, structures that have continuous rota- 
tional symmetry about an axis (Coo) are not chiral. The 
simplest nonuniaxial structures have second-rank mass- 
moment tensors with three inequivalent principal axes 
and are biaxial. Fig. 2 shows model linear, biaxial pla- 
nar, and chiral molecules. The twisted "H" in Fig. 2c is 
both biaxial and chiral. 




(a) (b) 




(c) (d) 

FIG. 4. Examples of chiral structures created from achiral 
ones: (a) left, an achiral molecule in the shape of a tetrahe- 
dron with four equal masses at its vertices; right, a similar 
chiral molecule with four unequal masses at its vertices, (b) 
left, an achiral planar sheet; right, a helix formed by twisting 
a sheet about an cylinder, (c) left, an achiral propeller with 
all blades perpendicular to the hexagonal core; right, a chiral 
propeller with all blades rotated away from the normal to the 
hexagonal plane, (d) left and middle, cylinder with achiral 
decorations; right, chiral cylinder with helical decorations. 

In spite of its practical importance, there is no uni- 
versal quantitative description of molecular chirality nor 
is there an accepted procedure for identifying the chi- 
ral part of an intermolecular interaction. As a result, 
only recently has real progress been made in addressing 
straightforward and apparently simple issues such as the 
relation between the cholesteric pitch and molecular ge- 
ometry. Ideally one would like to introduce a parameter 
that measures the chiral strength or degree of chirality of 
a given molecule. A non-vanishing value of this parame- 
ter, which we will refer to as a chiral strength parameter 
or simply a chiral parameter, would distinguish a chiral 
molecule from an achiral one, just as the dipole moment 
distinguishes a polar molecule from a non-polar one. In 
addition, such a chiral parameter would play a crucial 



role in determining macroscopic chiral properties, such as 
the optical rotary power, the wavevector for cholesteric 
ordering, and other macroscopic chiral indices. Unfortu- 
nately, there appears to be no such simple description of 
chirality and chiral parameters. As we will show, just as 
a charge distribution can be described by an infinite hier- 
archy of multipole moments, so a chiral molecule can be 
described by an infinite hierarchy of chiral parameters. If 
any one of these parameters is nonzero, then the molecule 
is chiral. Also, since different macroscopic properties will, 
in general, depend on different microscopic chiral param- 
eters, we do not expect strong correlations between the 
various macroscopic manifestations of chirality. 




(b) 

FIG. 5. This figure shows how steric interactions between 
two chiral molecules produce a net relative rotation of their 
long axes. The "barber pole" stripes on the model cylin- 
drical molecules represent protrusions such as are found on a 
screw. If the protrusions on one molecule (the shaded regions) 
are forced to fit "hand-in-glove" into the grooves on another 
(the unshaded regions) then, as we show on the right, when 
one molecule is placed upon another, the protrusions on one 
molecule align with the grooves on the other. The result is 
that the long axes of the two molecules acquire a relative 
twist determined by the pitch of the stripes. In (a) we show 
molecules with a tight right-handed pitch as determined by 
the geometrical right-hand rule. The relative twist of two 
neighboring molecules is right-handed according to the right- 
hand rule, (b) shows molecules with a weak right-handed 
pitch. The relative twist of neighboring molecules is now left- 
handed. These examples show that the "handedness" of in- 
dividual molecules does not determine the handedness of the 
collective structure. 
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Even if one were equipped with a complete under- 
standing of the nature of interactions between chiral 
molecules, the calculation of macroscopic parameters like 
the cholesteric pitch is not completely straightforward. 
An argument due to Straley (1974) and illustrated in 
Fig. 5 makes it clear why molecular chirality causes 
macroscopic rotation: when two screw-like molecules are 
brought close together, their grooves interlock to pro- 
duce a finite rotation angle A9 between long molecular 
axes. A simple estimate of the pitch based on this pic- 
ture is P rts (2tt/A9)1 where I is the molecular diame- 
ter. Taking a rough estimate of AO w 10° — 20° and a 
molecular size of roughly 1 nm, one finds that P w 10 
nm, two or three orders of magnitude smaller than typ- 
ical pitches. Indeed some chiral systems (Fraden, 1995) 
are labeled nematic rather than cholesteric, presumably 
because their pitches are too long to be measured exper- 
imentally. Thus a quantitatively correct calculation of 
the cholesteric pitch cannot be based solely on molecular 
parameters and presents a challenge to theorists. As we 
demonstrated previously (Harris, Kamien, and Luben- 
sky, 1997), for central-force or steric models, go vanishes 
(infinite pitch) unless biaxial correlations between the 
orientations of adjacent molecules (such as are illustrated 
in Fig. 1) are taken into account. Quantum interactions, 
however, do not require such correlations and hence can 
give a nonzero value of go even within mean field theory. 

The purpose of this paper is to present recent progress 
both in quantitatively characterizing molecular chirality 
and in calculating the cholesteric pitch from microscopic 
interactions. In Section II we start by making some naive 
qualitative comments about the nature of chiral symme- 
try. The central idea is that an achiral object has higher 
symmetry than a chiral one. We will develop a systematic 
procedure for generating a countably infinite set of chiral 
molecular parameters that all vanish when the molecule 
is achiral. Next, in Section III, to illustrate our chiral 
parameters we consider a topological "rubber glove" , a 
chiral structure that can be converted to its mirror im- 
age via distortions through a continuum of intermediate 
states all of which are chiral. This demonstration shows 
clearly that handedness is not an absolute concept, but 
depends on the property under consideration. In Section 
IV we describe a calculation of the cholesteric pitch from 
a classical model of central forces between atoms. This 
calculation shows that biaxial correlations play a critical 
role in determining the pitch P. In fact, if these biaxial 
correlations do not exist, each molecule rotates freely and 
appears, on average, uniaxial and thus achiral. In accord 
with our discussions of chiral parameters in Section II, we 
expect that other macroscopic chiral response functions, 
such as the rotatory power, will depend on other chiral 
structure parameters. Indeed, it is likely that an under- 
standing of many such indices of molecular chirality is 
required to interpret the dramatic frequency dependence 



of these susceptibilities. 

II. CHIRAL PARAMETERS 

As mentioned in the introduction, one expects the chi- 
ral interaction between molecules to involve parameters 
characterizing the chiral strengths of the molecules. How- 
ever, there is no obvious precise quantitative formulation 
of parameters that characterize the degree to which a 
given molecule is chiral. Only a handful of chiral strength 
parameters have been proposed. For instance, Osipov et 
al. (1994) developed a measure of molecular chirality by 
considering the symmetry of response functions describ- 
ing the electromagnetic behavior of chiral molecules. We 
also introduced a chiral strength parameter in a previous 
calculation (Harris, et. al, 1997) of the chiral wavevector 
go- We will review this calculation of go m Section IV. 
Both of the above chiral parameters are non-local (in a 
sense to be made more precise later) - a feature which we 
shall see is generic. In the following we give a more sys- 
tematic discussion of the structure of such chiral strength 
parameters. 

A. Chirality is the Absence of Symmetry 

A preliminary remark is that "chiral symmetry" is ac- 
tually an absence of symmetry, i.e., the absence of sym- 
metry under improper rotations. Thus a chiral object 
has lower symmetry than an achiral object which is in- 
variant under chiral operations. This situation contrasts 
with "spherical symmetry" which implies the existence, 
rather than the nonexistence, of symmetry elements. Ac- 
cordingly, it is instructive to compare the way spherical 
symmetry is destroyed when a sphere is distorted to the 
way achiral symmetry is destroyed when an achiral object 
is chirally distorted. To start, we consider distortions of 
a sphere centered at the origin. Initially, the sphere of 
radius ro is a surface described by 

r = r , (2) 

where r is the radial coordinate from the origin. When 
the sphere is distorted, the radial coordinate of its sur- 
face will depend on the usual angles 9 and 4> and can be 
expanded in spherical harmonics: 

oo m—n 

r = r + ^2 ^2 a nm Y nm {9, (f>) , (3) 

n— 1 m=—n 

where a nm = a* _ TO (-l) m . 

Usually one characterizes distortions by the values of 
the 

o-nm 's for the smallest value of n for which one of these 
is nonzero. The first-rank tensor a\ m describes transla- 
tions of the sphere, which do not alter the symmetry 
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and which we ignore. Thus, the lowest-order distortions 
are characterized by ci2 m , which in a Cartesian repre- 
sentation is a symmetric, traceless, second-rank tensor. 
In general, a complete specification of the shape of an 
asphcrical surface requires the values of the infinite set 
of a nm . Since the a nm mix with each other under ro- 
tation, it is desirable to construct rotationally invariant 
measures of asphericity. A useful class are the quantities 

<7 n =J2a 2 nm , (4) 

m 

which provide rotationally invariant characterizations of 
the magnitude of the asphericity associated with nth- 
rank tensor distortions. It is entirely possible for 02 to 
vanish while higher-order cr„ do not. In this case the dis- 
tortion is characterized by the lowest-order, nonvanishing 
value of a n . 

With the above discussion in mind, we consider chi- 
rality. An object can be described by the multipole mo- 
ments of its density, p(r, 9, <p), namely 

PimN = J dxp{r,8,<l>)T»Y lm {6,<l>) , (5) 

where p* lmN = {-l) m p~i- m ,N since p(r,9,<p) is real. The 
moments p~i m N for a given N define a tensor parame- 
ter that transforms under a (21 + l)-dimensional repre- 
sentation of the rotation group. The alternative Carte- 
sian representation in terms of symmetric, traceless ten- 
sors pYn °f ran k I is use d extensively in treatments of 
liquid-crystalline order, and we will employ them when 
appropriate. For a molecule consisting of point atoms, 
the density consists of a sum of Dirac-delta functions lo- 
cating each atom. This provides a natural framework to 
study interactions between molecules, in which connec- 
tion the central quantities are p~i m N- The question we 
wish to address here is how these moments, or appropri- 
ate functions of them, characterize chirality. We start by 
discussing the analogs of the parameters a n in order to 
characterize the magnitude of chirality. 

B. Construction of Pseudoscalars 

Bearing in mind that chirality requires the absence 
of inversion symmetry, we propose to characterize the 
magnitude of chirality by an infinite sequence of pseu- 
doscalars. First, note that Lord Kelvin's definition may 
alternatively be stated as "an object is achiral if there 
exists a rotation Jl such that the object is invariant un- 
der the operation i~iS2, where S2 is spatial inversion." 
Since spatial inversion is a mirror operation followed by 
a 7r rotation about an axis perpendicular to the plane 
of the mirror, this definition of chirality is equivalent to 
Lord Kelvin's. A scalar is invariant under both rotations 



and inversion while a pseudoscalar is only invariant un- 
der rotations - it changes sign under inversion. Thus any 
pseudoscalar parameter ij) n that we construct from the 
multipole moments of the density will necessarily vanish 
when the molecule is achiral. Furthermore, the degree of 
chirality can be characterized by the magnitudes of the 
set of ifj n just as the degree of asphericity was character- 
ized by the various a n . 

The construction of pseudoscalars can be done system- 
atically by considering the representation theory of the 
three-dimensional rotation group 0(3). This procedure 
amounts to nothing more than considering the quantum- 
mechanical addition of angular momentum. To each rep- 
resentation we will attach its transformation properties 
under inversion (i.e., parity). Pseudoscalars will trans- 
form as 1-dimcnsional representations with odd parity. 
Denoting the <i-dimensional representation with parity p 
as cP, we note that the rank-Z representations generated 
by the multipole expansion are 

l+,3-,5+,...,(21 + l) ( - )! ,.... (6) 

To construct a pseudoscalar we must form tensor prod- 
ucts of different representations. While the resulting rep- 
resentations obey the rules for addition of angular mo- 
mentum, the parity of the new representation is simply 
the product of the parities of the two representations. 
For instance, since two spin-1 states (with odd parity) 
can be combined to form a spin-2, spin-1 or spin-0 state 
(all with even parity), we have 

3- ®3" = 5 + ©3 + ©l+ . (7) 

This gives us our first representation that is not a mul- 
tipole representation: 3+ a pseudovector. Forming the 
triple tensor product: 

3 ® 3~®3~ = 

7~ ©5" ©5" ©3~ ©3~ ©3" 01", (8) 

we arrive at our first pseudoscalar 1 , which we recog- 
nize as the vector triple product A • (B x C). The above 
discussion suggests that any pseudoscalar must involve a 
product of at least three of the multipole moments pimN- 
It is, in fact, always the case that a 1~ can only be con- 
structed from the tensor product of d + © d~. Since the 
multipole moments do not include both d + and d , one 
must construct one of these via tensor products. Thus 
any pseudoscalar must involve a product of at least three 
of the multipole moments. This implies that a chiral pa- 
rameter for a given object can be expressed as an integral 
over at least three position vectors in that object, and in 
this sense, chirality is a nonlocal property. 

1. Nonpolar Molecules 

In this subsection we confine our attention to the case 
in which vector representations can be eliminated by 
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proper choice of the center of the molecule. Specifically, 
in this case we do not allow molecules to have a dipolc 
moment. (That case will be considered in the next sub- 
section.) To illustrate this theoretical discussion, let us 
look for the lowest order (in powers of components of 
r) pseudoscalar that can be constructed from the simple 
mass- weighted distance moments 

Pim = E \ r x\ l Ylm{0 X ^x) ( 9 ) 

of a homoatomic molecule, where the sum is over atoms \ 
in the molecule X. Throughout, we will label molecules 
with capital Roman letters and their constituent atoms 
by Greek letters. The center of mass of molecule A will 
be R A , and each atom a will be displaced from there by 
r Aa . For simplicity we focus only on pi m = pi m i- Our 
discussion could be embellished by considering pi m N for 
other values of N. If we measure the density relative 
to the center of mass, then p\ m = for all m. Thus the 
multipole expansion only provides us with d = 5 or larger 
dimensional representations. Since we have restricted 
ourselves to a single tensor for each (i-dimensional repre- 
sentation, 5+ ® 5+ will not contain any pseudotensors, 
and the lowest-order pseudotensor we can construct is 
contained in 5 + <8> 7~ . We could now try to construct 
a pseudoscalar by contracting the resulting tensors with 
5+, but, again because we are considering only moments 
of pi m , we would get zero just as we would get zero for 
the triple product A x B A = in the vector case. 
A nonzero pseudoscalar only results when 5 + <S> 7~ is 
contracted with a tensor different from the 5 + and the 
7~. Thus the lowest order pseudoscalar we seek is in 
5 + ® 7~ ® 9 + . In terms of spherical harmonics we set 

■00 oc E C(234; mn)p 2m P3np4, m +n i ( 10 ) 

mn 

where C(234;mn) are the appropriate Clebsch-Gordon 
coefficients. It is convenient to choose the normalization 
so that in the Cartesian representation, this is 

%=pl 3 pf m e ikp pi lmp • (11) 

Note the presence of the antisymmetric symbol e^fe in 
this expression. It is required to produce a scalar from 
two even-ranked and one odd-ranked tensor. 

We can calculate ipo for the "twisted H" molecule, Mi, 
shown in Fig. 2, with four identical atoms at 

Mi = {(a, 6, c), (a, -b, -c), (-a, b, -c), (-a, -b, c)} . 

(12) 

This molecule is chiral if abc ^ and if \a\ ^ \b\ ^ \c\ ^ 
\a\. Since 0o is a rotational invariant, it may be evaluated 
in any convenient coordinate basis. We find 

0o = K (a 2 - b 2 )(b 2 - c 2 ){c 2 - a 2 )abc , (13) 



where Kq is a numerical constant. Note that tpo does 
indeed vanish when the parameters assume the special 
values for which the molecule has the higher achiral sym- 
metry. Moreover, since the mirror image of Mi may be 
obtained by exchanging any two of a, b and c or by re- 
versing any one of their signs, we see that under inver- 
sion ipQ — ► —-00 , an d it is indeed a pseudoscalar. This 
ninth-order multinomial is the lowest-order expression 
constructed from pi m , which must vanish for an achiral 
object. However, just as in the discussion of asphcric- 
ity, it is possible to consider a class of chiral molecules 
for which tp vanishes but which requires an even higher- 
order multinomial to characterize its chirality. Consider a 
twelve atom molecule, M 2 , obtained by taking the atoms 
as in the "twisted H" together with the eight atoms ob- 
tained by cyclic permutation, so that identical atoms are 
now at 

M 2 = ^(a,b,c),(a,-b,-c),(-a,b,-c),(-a,-b,c), 

(b, c, a), (-6, -c, a), (6, -c, -a), (-6, c, -a), (14) 
(c, a, b), (-c, a, -6), (-c, -a, 6), (c, -a, -6) j. 

To show that tpo vanishes for M 2 , it suffices to verify that 
pim = for all to. (This result is most easily verified 
using the Cartesian representation for p% ■) However, ex- 
cept for the special values of the parameters (i.e., a = 0, 
b = 0, c = 0, \a\ = \b\, \b\ = \c\, or |c| = |a|), this object is 
clearly still chiral since it is the union of three identical 
chiral objects. To describe its chirality the lowest order 
pseudoscalar constructed from the moments pi mi is 

i>i=^ C(346; ^)p^piuPl^ +v , (15) 

which we evaluate to be 

0! = K ia bc(a 2 - b 2 )(b 2 - c 2 )(c 2 - a 2 ) 

x (a 4 + 6 4 + c 4 - 4a 2 6 2 - 46 2 c 2 - 4c 2 a 2 ), (16) 

where K\ is a numerical constant. As was the case for tpo, 
this chiral strength parameter vanishes when the "twisted 
H" is made to be achiral. 

It is clear that we can construct an infinite sequence 
of chiral parameters from the pi m that vanish for achiral 
objects. For example, when J, K, and L are all different 
integers whose sum is odd, each of the quantities 

IpJKL = C(JKL; mn)p Jm p K nP*L, m +n ( 17 ) 

mn 

is a chiral parameter, which is nonzero only for chiral 
molecules. 

If we consider different tensors of a given rank, we 
can construct other sets of chiral parameters not encom- 
passed by ipJKL- For example, if there are two distinct 
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second-rank tensors ^ 2 and t 2 , then we can construct 
the chiral parameter 



V2=l2 T 2 e vkP 3 

cc^2 C(223; mn)j 2m T2nP3, m+n - 



(18) 



When j 2 — r 2 , 4>2 vanishes because eijk is antisymmet- 
ric in all indices. The tensors 7 2 J could, for example, 
be constructed from p 2m arid p 2m N for N =/= I. Alterna- 
tively, two different tensors can be constructed from the 
second-rank mass-moment tensor 



ij 
P2 



E 

xex 



r Yx 



(19) 



To do this, we express p 2 in the basis of its orthonormal 
principal axes emblazened on the molecule ei, e 2 , and 
e 3 , where e 3 = ei x e 2 is associated with the largest 
magnitude eigenvalue of p 2 . Then we decompose p 2 
into its uniaxial (Q y ) and biaxial (_B y ) components as 

p 2 j = SQ ij + B (e\e{ - e 2 e J 2 ) = SQ ij + B ij , (20) 

where 



I \2 1 2 

(r x -e 3 ) --t x 



^E [( r x ' e i) 2 - Ox - e 2Y 



and 



(e 3 e 3 8 



(21a) 
(21b) 

(21c) 



Note that B % i , the biaxial part of p 2 , vanishes when the 

molecule is uniaxial. Setting ^2 — Q*- 7 ', t 2 = B^ , and 
p ijk = ^jfe^ wc obtain 

^ 2 = Q il Bi m e ijk S klm (22) 

as a chiral strength parameter, which plays a role in our 
calculation of the cholesteric pitch go to be presented in 
Section IV. For the "twisted H" molecule, we have 

i) 2 = C 2 abc [\a 2 + c 2 - 2b 2 \ + a 2 - c 2 ] 

for a 2 < b 2 < c 2 , (23) 

where C 2 is a numerical constant. Expressions for ip 2 
in regimes other than a 2 < b 2 < c 2 can be obtained by 
suitably permuting variables. One may verify that when 
Eq. (23) is valid, ip 2 does vanish when the molecule is 
achiral, i.e. when a 2 = b 2 < c 2 or a 2 < b 2 = c 2 

One can construct different second-rank tensors in 
other ways. For instance, in some phenomenological in- 
tcrmolecular potentials, the strength of the dispersion 



interaction between a pair of atoms is estimated to scale 
with the product of their atomic polarizabilities. In that 
case, a polarizability-weighted second-distance moment 
is generated by the multipole expansion of the inter- 
molecular potential. In general, a tensor with any weight- 
ing that is distinct from the mass weighting can play the 
role of the additional second-rank tensor needed to char- 
acterize chirality. Such moments would have the form 



Plr, 



= Yl W X Y lrn(0x,<t>x) ' 
X6x 



(24) 



where w x is a weight factor, which can differ from the 
factor r x in Eq. (9). For a molecule with p atoms, we 
can obviously have up to p linearly independent second- 
rank tensors. As we mentioned, moments similar to these 
have been used in the study of optical properties of chiral 
systems by Osipov et. al (1995). 



2. Ferroelectric Liquid Crystals 

There are cases in which one may invoke the vector 
representation, even in liquid crystalline systems. A 
particularly interesting case is that of ferroelectric liq- 
uid crystals (Meyer, et al., 1975). These phases are 
composed of mesogens that have an electric dipole mo- 
ment. Recall that the smectic-C liquid crystalline phase 
is a one-dimensional layered structure with layer normal 
N. In each layer the nematic orientation n is not par- 
allel to N. Thus, one can construct the pseudovector 
A = (n • N)(n x N). It is clear that under parity A docs 
not change sign and that both the nematic (n — > — n) 
and layer normal (N — > — N) inversion symmetries are 
preserved. However, if the molecules are chiral then, as 
we have seen, a nonvanishing pscudoscalar ip may be con- 
structed. In this case P = ipA is a true vector and can 
set an unambiguous alignment direction (perpendicular 
to n and N) for the molecular dipole moments. 

In general, when electrostatic interactions are taken 
into account, both signs of charge are present and there 
is a non-zero dipole moment that no change of origin 
can eliminate. In this case, we can again construct two 
different second-rank tensors t 2 and ~f 2 , and then 

4>f - Yl C ( 122; mn )PlmT 2 nl 2 . m+n OC (^4^2^31, (25) 
mn 

is a pseudoscalar which includes the (dipole) charge mo- 
ment pi m . 



Pi 



E 

xex 



(26) 



In fact one could construct a pseudoscalar from a vector 
triple product of three noncoplanar vectors obtained by 
introducing three different weight factors into the sum in 
Eq. (26). 
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III. "RUBBER GLOVE" MOLECULES 

We have argued that a quantitative characterization 
of chirality does not rest on one parameter, but rather 
on an infinite hierarchy of chiral moments. However, one 
has a natural tendency to associate a specific handed- 
ness to a given chiral molecule. We will show that even 
the "handedness" of a molecule depends on the chiral 
property under consideration. This is, in fact, familiar 
from circular dichroism measurements: the difference in 
attenuation of left- versus right-circularly polarized light 
changes sign as a function of its wavelength. Thus the 
handedness of an object is really in the eye of the be- 
holder. 

TABLE I. Atoms and weightings W and W' for the chiral 



parameters ipo and ip' of the "twisted H" molecule, Mi . 



Position (x) 


W( X ) 


W'{x) 


(a, b, c) 


1 


1 + M 


(a, -b, -c) 


1 


1 + p 


(—a, b, — c) 


1 




(-a, -b, c) 


1 





To illustrate this idea, we consider a process in which 
a chiral molecule is continuously deformed into its enan- 
tiomer or mirror-image molecule. For the "twisted H" 
molecule, we could do this by continuously varying the 
parameter a until its final value becomes the negative of 
its initial value. Obviously, when a passes through zero 
the molecule is achiral and one might be tempted to say 
that the plane a = in parameter space separates re- 
gions of opposite handedness. However, as we will show 
by example, it is possible to continuously deform a chi- 
ral molecule into its enantiomer without passing through 
an achiral configuration. (Here "deformation" is used in 
its most general sense in which not only the coordinates, 
but also the mass and other properties of atoms are var- 
ied.) The existence of such a continuous deformation is 
incompatible with the existence of an intrinsic definition 
of right or left handedness. It is also clear that any sin- 
gle measure of chirality will pass through zero at some 
point in the process of deforming a molecule into its enan- 
tiomer. However, a molecule is achiral only if all of its 
chiral moments are zero; the vanishing of a single chiral 
moment alone does not make a molecule achiral. We will 
illustrate explicitly that there exist paths of deformations 
between enantiomers along which there is no point where 
all chiral moments vanish. Nevertheless, along this path 
every chiral measure must and does vanish at some point. 
A molecule that can be deformed in this way is known as 
a topological rubber glove in analogy with a real rubber 
glove - it can be inverted one finger at a time thus always 
remaining chiral (Walba, et al, 1995). In the context of 



our discussion we would interpret this by saying that the 
eye automatically measures many indices of chirality, and 
as each finger is inverted some indices may pass through 
zero to change sign until finally all indices have changed 
sign. 

We can see this explicitly by keeping track of more 
than one chiral parameter as the "twisted H" is inverted 
continuously into its enantiomer. The two chiral param- 
eters we will monitor are ipo and i()' . ipo is defined in Eq. 
(10) and given explicitly in Eq. (13). ip' is also defined as 
in Eq.(10), except that now pi m is replaced by pi m \W'\, 
where 

Pi m [W] = J2 W(x)\r x -r°[W]\ l Y lm (6 x ,ct> x ) , (27) 

where W{\) is a weighting function associated with 
the atom x an d r°[TT] is the IT- weighted center of 
the molecule, chosen so that pi m [lT] vanishes. Various 
weighting functions are shown in Table I. Until now we 
have considered molecules composed of identical atoms, 
i.e. W(x) = I- Of course, the molecule need not have 
identical atoms, and, therefore, not all properties of the 
atoms need be the same. For instance, if all the atoms are 
weighted equally, p/ m [IT] would correspond to a purely 
geometric moment. However, if we were to weight the 
atoms by their polarizabilities, then the moments would 
be different. We have defined Vo to reflect geometric 
properties and tl>' to reflect others. We could have in- 
stead introduced other weight functions which reflect the 
valence, electronegativities, etc., of the atoms. Since dif- 
ferent properties are not perfectly correlated, they may 
require different weight functions. To construct the con- 
tinuous deformation between enantiomers, we only in- 
voke "twisted H" molecules that have atoms with polar- 
izabilities 1 + jU on sites #1 and #2 and 1 — p, on sites 
#3 and #4. The point here is that the molecule is chiral 
if either ?/>o or ip' is nonzero. 

In the calculation of ip' all displacements r x are eval- 
uated relative to the "center of polarizability" of the 
molecule, so that p' lm = for all m. We find that 

V>o = K (l -p 2 ) 2 abc{b 2 -c 2 ) 

x[(c 2 -a 2 )(a 2 -b 2 )-lp 2 a 4 } . (28) 

When the molecule is tetrahedral (a = b = c), it is truly 
achiral if any two of its atoms are identical. Indeed in this 
case tpo and ip' (as well as all other chiral parameters) 
vanish. 

We will now consider a process in which the molecule is 
distorted in the parameter space (a, b, p) from the initial 
configuration A of Fig. 6 into its enantiomer, E, while 
remaining chiral along the entire path of deformation. 
Initially p = and < a < b < c. 
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b(0) 



a(0) 



a(0) 
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FIG. 6. Path I (AFE) and path II (ABCDE) between chiral 
enantiomers (mirror images) A and E for a molecule consisting 
of 4 atoms at the vertices of a twisted H for the parameters 
used in Fig. 7. Below each configuration we give (a, a'), 
where a and a' are respectively the signs or zero value of ip 
and ip' . In configurations C and F, a = b. Path AFE passes 
through the achiral configuration F. Path ABCDE passes 
through chiral configurations only. Configuration C would 
have a mirror plane (as does F) if the masses were all equal. 

Note that any pseudoscalar must change sign under 
reflection and therefore must pass through somewhere 
along the path between enantiomers. We will consider 
two paths between the molecule and its mirror image 
described in Fig. 6: the first will be a path AFE through 
an achiral point, the second a path ABCDE that goes 
only through chiral states, along which ifj and ip' Q never 
simultaneously vanish. Our deformation will rearrange 
the molecule into its mirror image under the operation 
(x,y,z) —> (y,x,z) which takes (a, 6, c) into (6, a, c). We 
may parameterize AFE by 



a(t) = a(0) + [6(0) - a(0)]t 
b(t) = 6(0) + [o(0) - 6(0)]t 
li(t) = , 



(29a) 
(29b) 
(29c) 



where t = corresponds to the initial configuration A 
and t = 1 to the enantiomer E. For concreteness, we 
have used the values a(0) = 0.99, 6(0) = 1.01 and c = 1.2. 
Note that t = 1/2 corresponds to the point F, which is 
achiral, at which it is easy to see that Vo = V'o = 
since a(l/2) = 6(1/2) and tt = 0. To pass between 
enantiomers without passing through an achiral configu- 
ration we will follow the path ABCDE along which the 
mass parameter /u, does not remain fixed at zero. Note 
that this path avoids the line, a = b and n = 0, along 



which the molecule is achiral. Over the first section of 
the path, AB, we change the masses of the atoms by 
changing \i from its initial zero value to a suitable value 
/i (po = 0.15). Over the second section, BCD the mass 
parameter is held fixed, so that jj, = no, but a and 6 are 
varied as in Eq. (29a, b), so at D ^ = /i , a{to) = b 
and b{tr)) = a. Finally along DE, a(t) and b(t) remain 
constant but jj, is changed from /io back to zero. For 
the parameter values we have chosen, ip' changes sign on 
this part of the path. Since configuration E is the mirror 
image of configuration A, the values of all chiral param- 
eters, including ip and tj)' have changed sign. But all 
states in the path of deformation are chiral and nowhere 
on this path do both V'o and ip' simultaneously vanish! 
The chiral measures V'o an d Vo f° r tne path ABCDE 
are plotted in Fig. 7. 







ST o 








1 \ 




1 

B 




v D \ 


h- 

E 



FIG. 7. Plot of the chiral measure ipo and ip'o along the path 
ABCDE. Note that both parameters pass through zero on 
this path, but they do not pass through zero at the same place. 
The parameters used to obtain this plot are a(0) = 0.99, 
6(0) = 1.01, c = 1.2, and fi = 0.15. 

We note that in a molecule with more than four atoms, 
our artificial deformation of the masses can be replaced 
by the additional degrees of freedom provided by the 
other atoms. 

What shall we conclude from this example? Since we 
can continuously deform a molecule into its enantiomer 
via only chiral states, there is no general, unambiguous 
characterization of handedness. 



IV. PREDICTION OF THE CHOLESTERIC PITCH 

We now turn to the calculation of the cholcsteric 
wavenumber qo in terms of the microscopic interactions 
between molecules. We will argue that previous clas- 
sical analyses of this problem have missed an essential 
feature of chiral interactions - the necessity of biaxial 
correlations between the molecules. Our result will in- 
volve the chiral parameter ip2 defined in Eq. (22). In ac- 
cordance with experimental observations we will assume 
that go a -C 1, where a is a typical intermolecular sepa- 
ration. This means that the cholcsteric can be treated 
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locally as a nematic even though it is twisting on longer 
length scales. Our aim, then, is to obtain a formula for 
(jo in terms of correlation functions evaluated in the ne- 
matic limit, i.e., when all chiral interactions have been 
turned off. This type of result is particularly desirable for 
numerical simulations, since it only requires nematic cor- 
relation functions and hence simulations of the nematic 
state rather than the cholesteric state. 

In this spirit we start by considering the long- 
wavelength properties of systems that are locally ne- 
matic with fluctuations described by the phenomenolog- 
ical Frank free energy (Frank, 1958), which depends on 
the director n(x). When surface terms are neglected, 
which only come into play when topological defects or 
internal surfaces are present (Kleman 1983), the Frank 
free energy is 

F=^Jd 3 x (V • n) 2 + K 2 (n • V x n) 2 

+K 3 [n x (V x n)] 2 + 2/in ■ V x nj. (30) 

This effective free energy can describe most of the phe- 
nomena of nematic and cholesteric liquid crystals. The 
parameter h is the generalized thermodynamic force that 
determines the pitch. We shall refer to h as the torque 
field since it is proportional to the microscopic inter- 
molecular torques as we shall show. To see how h de- 
termines the pitch, we consider a helical configuration of 
n(x) as in Eq. (1). Then, the twist 

-n-Vxn = <? (31) 

is spatially uniform, V • n = 0, and n x (V x n) = 0. 
The Frank free energy is 

F = T^K 2 q 2 -h q y (32) 

where T is the volume of the system. This energy is 
minimized when 

q = qo = h/K 2 . (33) 

Note that ft. is a pseudoscalar. Under spatial inver- 
sion, h and, therefore, also the "twist" change sign: 
n • V x n — > — n-Vxn since it is linear in spatial gra- 
dients. To obtain a nonzero value of h it is, therefore, 
necessary that the system not be invariant under spatial 
inversion. Furthermore, we see that to evaluate qo we 
need to calculate both h and K 2 . 

The Frank elastic constants K\, K 2 , and K 3 can be 
estimated by dimensional analysis using only excluded 
volume (cntropic) interactions by dividing an energy 
scale by a length scale. Taking the energy scale to be 
k B T ~4x 10 _14 erg and the length scale to be a molec- 
ular size ~ 1 nm, we find that Ki ~ 4 x 10~ 7 dync, or 



about 1 fi dyne. This is approximately correct for typi- 
cal liquid crystals that exist at room temperature. As we 
discussed in the introduction, dimensional analysis does 
not predict qo correctly: qQ 1 is typically on the order of 
or larger than hundreds of nanometers and does not cor- 
respond to any natural length scale in the problem. The 
challenge, therefore, is to calculate h and to determine 
why its magnitude does not correspond to what dimen- 
sional analysis suggests. 

A major objective is to establish an approach that 
in principle will provide a rigorous calculation of the 
cholesteric pitch (or cquivalently of h) in the limit when 
the pitch is very long, i.e., when q$a — > 0. The start- 
ing point of any microscopic calculation of h must be the 
intermolecular potentials. Liquid crystal mesogens are 
notoriously complex, containing hundreds to hundreds 
of thousands of individual atoms, and they have corre- 
spondingly complex interactions. A reasonable approach, 
and the one we will pursue here, to construct the de- 
sired potentials is to model each mesogen as a collection 
of connected spherically symmetric atoms that interact 
via pairwise central-force potentials with atoms on other 
mesogens. The interatomic potential consists of a long- 
range van der Waals part and a short-range repulsive 
part arising mostly from the Pauli exclusion principle. 
Fluid physics is dominated by the short-range repulsive 
part, and its is often useful to replace the full interatomic 
potential by a simple hard-core potential with no attrac- 
tive part. An intermolecular potential constructed in this 
way includes steric interactions that force two chiral ob- 
jects like screws or ridged "barber poles" (Fig. 5) to twist 
relative to each other when in contact. 

There are contributions to the effective intermolecu- 
lar potential that cannot be expressed as a superposi- 
tion of interatomic central-force potentials. The simplest 
such contribution is a chiral dispersion potential - an 
anisotropic generalization of the van der Waals poten- 
tial. Dispersion forces arise from the Coulomb potential 
between all pairs of electronic and nuclear charges and 
quantum fluctuations of the electronic states. The van 
der Waals potential is produced by the interaction of fluc- 
tuating electric dipoles on different atoms. If molecules 
are chiral, a fluctuating dipole on one molecule can inter- 
act with a fluctuating quadrupole on another to produce 
an effective chiral dispersion force (van der Meer et al, 
1976; Kats, 1978; Issaenko et al, 1998). It is difficult 
at the moment to obtain first-principles estimates of the 
strength of chiral dispersion forces. 

In what follows, we confine our attention to pairwise 
central-force interactions. Once the intermolecular po- 
tentials have been chosen, the next step is to devise a 
scheme to compute h. Since q a <C 1, one may assume 
that h is small and calculate all quantities to lowest or- 
der in h, or equivalently to lowest order in q a. The 
cholesteric twist induces biaxial contributions to the nc- 
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matic order parameter or order (qoa) 2 (Priest and Luben- 
sky, 1974). Thus, to lowest order in q^a, biaxiality can 
be ignored, and the cholesteric can be treated as though 
it were locally uniaxial. Mean-field theory is a natural 
first calculational approach to pursue (Schroder, 1979; 
Evans, 1992; Pclcovits, 1996; Moro, et. al, 1996). In 
the locally uniaxial limit appropriate to most cholester- 
ics, mean-field theory will always predict h = (Salem 
el al. 1987) when central-force potentials between atoms 
are assumed. This result is easy to understand: Mean- 
field calculations seek the best self-consistently deter- 
mined distribution function for a single mesogen. In a 
uniaxial system, this distribution function will be uni- 
axial and produce only uniaxial average mass moments. 
Since there are no uniaxial structures that are chiral, any 
manifestation of chirality is washed out, there will be no 
potential favoring relative twist of neighboring molecules, 
and go will be zero. Thus, a more powerful approach than 
mean-field theory is needed to calculate qo in the majority 
of cholesterics that are nearly uniaxial. In the less com- 
mon case that would arise when chirality is introduced in 
a biaxial nematic, the cholesteric is locally biaxial, and 
mean-field theory will produce a nonvanishing value of 
qo- 

The failure of mean-field theory can be traced to 
its neglect of biaxial correlations between neighboring 
molecules. A first principles theory developed by the 
authors (Harris et al., 1997) provides a rigorous method, 
not limited to mean-field theory, for calculating h. Its 
principal result is that, under certain approximations, 
h is proportional to a measure ip of molecular chirality 
times the spatial integral of a biaxial correlation function 
- a function which is strictly zero in mean-field theory. 
Thus, h is small and deviates from expectations based on 
dimensional analysis both because ip can be small and 
because biaxial correlations may be very short-ranged. 
Here we outline some important features of this theory. 
It begins with a rigorous expression of h, which can be 
obtained from Eq. (32): 



h = - 



1 dF 
T ~dq~ 



(34) 



9=0 



In the Appendix we show that within certain simplifying 
conditions this formulation leads to the result for the 
cholesteric wave vector, 



90 = ~4^T\ 



(35) 



where t ba is the torque exerted on molecule B by 
molecule A: 

t ba = e vkr J B pd k V (R + r Bfj - r Aa ) , (36) 

13a 

and (•) denotes thermodynamic averaging. It is no sur- 
prise that the intermolccular torques are the origin of 



the cholesteric structure. Indeed, if we had considered 
two planes of molecules a distance L apart along an axis 
perpendicular to the nematic director, then the change 
in angle between them would be 9 = qL and thus the 
expression for h in Eq. (34) would become 



— — 

T ~~d0 



(37) 



0=0 



which is simply the torque per unit area. Moreover, Eq. 
(35) provides a rigorous small-go expression for go in a 
fully aligned nematic in terms of quantities to be evalu- 
ated in the nematic limit, i.e., when molecular chirality 
is turned off. It is interesting to observe that this re- 
sult does not involve simply the torque r A = ^2 B t ba on 
molecule A. In the nematic phase the average torque on 
molecules in the interior of the sample is zero. Whether 
or not the nematic is locally unstable relative to states 
with nonzero q depends on the boundary conditions (Har- 
ris et al., 1999). Accordingly, Eq. (35) involves what 
we call the "projected torque on molecule A," namely 
^2 B H± ■ t ba . Finally, the appearance of the antisym- 
metric tensor in Eq. (36) guarantees that go is a pseu- 
doscalar and hence must vanish for a system in which all 
molecules are achiral. 

We now discuss some of the implications of the result 
in Eq. (35). For that purpose we consider a number of 
approximations that lead to a simple, yet nontrivial case. 
First, the molecules were assumed to have their long axes 
perfectly aligned along the director, i.e. their principal 
axes vectors e 3 are parallel to n. Although strictly speak- 
ing this limit is not realized in real systems, it does enable 
us to see some simple consequences of our formalism. Sec- 
ond, we neglect correlations between density fluctuations 
and orientational fluctuations. Finally, we will invoke an 
expansion in powers of r/R. Elsewhere (Harris et al., 
1999) we give a less restrictive discussion. As our prior 
discussion indicates, we must be sure to take the uniaxial 
average of the molecular orientations. On doing so, we 
found (Harris et al., 1997) that the first non-zero term 
upon averaging was fifth order in powers of r. The lowest 
order expression for go depends on the nematic alignment 
tensor Q l i , the biaxial tensor , and the third-rank ten- 
sor S ljk introduced in Eq. 21. We found 



qo = — 



E BA e ijk Q ip ((Bi l Sl pl 



B A l S k f 



8K 2 T 



(38) 



where K(R) is a sum of products of |R| and derivatives 
of the interaction potential V(R). 

The sum in Eq. (38) is averaged over molecular ori- 
entations and locations. The tensors B 1 ^ and S l ^ k de- 
pend on the orientation of the molecule. In fact, if the 
molecules spin independently about their long axes, this 
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average of B^ will vanish. Moreover, the components of 
S 1 J k that contribute to Eq. (38), , can be expressed 
in terms of B l £ so that, for identical molecules, Eq. (38) 
becomes 

= -fc ^ 2TTr( ^ 2) ■ (39) 

where ip 2 = S klm e ijk Q il B jm , as defined in Section II, 
is evaluated on a single molecule. Hence ip 2 is a pseu- 
doscalar parameter that is a measure of the chiral in- 
teraction between identical chiral molecules. It vanishes 
when the molecules are not chiral, and it provides a quan- 
titative index of chirality as would be measured through 
the cholesteric pitch. We emphasize, however, that other 
microscopic measures of chirality will in general involve 
other chiral parameters. Note that since the biaxial cor- 
relations can be negative at the intermolecular separa- 
tion, there is not even a correlation between the signs of 
ip 2 and the cholesteric pitch. 

The correlation function in Eq. (39) may be evaluated 
in the decoupling approximation where the molecular 
separation R is uncorrelated with the biaxial orientation. 
In this case the biaxial correlation function is simply the 
average 

((B«B«)) = (cos[2(0 A -^ B )]>, (40) 

where (see Fig. lc) 4> A is the angle between the biaxial 
axis of molecule A and the x-axis. 

In the above calculation of go, we assumed that all 
molecules in the cholesteric were chiral. In order for a 
chiral interaction producing relative twist between two 
molecules to exist, however, it is only necessary for one 
molecule to be chiral. The linearity of the our expres- 
sion for qo [Eq. (39)] in i\) 2 is a consequence of this fact. 
If both molecules had to be chiral, one might expect qo 
to be proportional to ip 2 - This is impossible, of course, 
because a pseudoscalar (go) cannot be proportional to 
the square of a pseudoscalar ("02 )• The expression for qo 
for a system composed of a mixture of chiral and achi- 
ral molecules is essentially the same as Eq. (39) except 
that A and B refer to different molecular species and 
if>2 is the chiral parameter of the chiral molecule. There 
are also chiral interactions between a chiral molecule and 
a strictly uniaxial molecule. These interactions lead to 
contributions to go that depend on correlations between 
the chiral parameter of the chiral molecule and the biax- 
ial anisotropy of the positional correlation of its uniaxial 
neighbor. Thus, when a nematic is doped with chiral 
molecules, a finite pitch must result, and we expect go to 
be proportional to the dopant concentration when it is 
small. 



V. DISCUSSION AND CONCLUSION 

In this paper we considered ways of characterizing and 
quantifying molecular chirality and of calculating the 
pitch wavenumber go, a macroscopic manifestation of chi- 
rality in cholesteric liquid crystals. We showed that there 
is not one, but an infinite number of chiral parameters 
that characterize a chiral object. Each chiral parameter 
is a pseudoscalar, whose construction, if it is obtained 
from mass or charge distributions, requires the contrac- 
tion of at least three mass- or charge-moment tensors. 
Chiral parameters for a given object can have varying 
magnitudes and even different signs. It is possible to 
pass continuously from a chiral object to its mirror im- 
age without ever passing through a state in which the 
object is achiral. In this process, each chiral parameter 
passes through zero, but at no point do all parameters 
pass simultaneously through zero. We showed in Section 
IV that the macroscopic pitch depends on both molecular 
chiral parameters and on molecular orientational corre- 
lations. Since these correlations vary with temperature, 
pH, pressure, etc., it is possible to change the magni- 
tude of the cholesteric pitch without changing molecules. 
This mechanism may be the explanation for the phenom- 
ena of twist inversion (Stegemeyer et al., 1989) in which 
the pitch changes continuously from right- to left-handed 
as a function of temperature. Moreover, orientational or- 
der (e.g. hexatic) can selectively enhance different inter- 
molecular correlations and thus change the importance 
of different chiral parameters for determining the macro- 
scopic cholesteric pitch. 

Finally, we emphasize the usefulness of the formal de- 
velopment of Section IV. It represents an important ad- 
vance in that in the limit of long pitch, it gives an ex- 
pression for the pitch in terms of quantities in the ne- 
matic system when chiral interactions have been turned 
off. Especially for simulations, this implies that it is not 
necessary to simulate a long pitch system. Instead one 
can simulate a homogeneous nematic in order to get cor- 
relation functions of the type appearing in Eq. (33). Fur- 
ther analysis in this direction may be needed to actually 
implement this idea. 
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GLOSSARY OF TERMS 

For the non-specialist, we include the following glos- 
sary of technical terms: 

1. cholesteric: a material in the cholesteric (or 
twisted nematic) phase. The director in this phase 
has a helical structure (depticted in Fig. 3) ob- 
tained obtained by twisting a nematic. 

2. director: unit vector n specifying the direction 
of average orientation of anisotropic molecules in a 
mesophase. 

3. enantiomer: a molecule with a given chemical 
formula can can exhibit many different geometri- 
cal structures called isomers. A chiral isomer is an 
enantiomer. A chiral isomer and its mirror image 
are an enantiomeric pair. 

4. Frank free energy: energy [Eq. (30)] associated 
with long- wavelength distortions of the director in a 
nematic. It is proportional to the square of spatial 
derivatives of n. 

5. mesogen: a molecule that forms a mesophase. 

6. mesophase: a phase with symmetry intermedi- 
ate between that of the most disordered isotropic, 
spatially homogeneous fluid phase and those of 
the most ordered period crystal phases. All 
liquid-crystal phases except for those, such as 
the cholesteric blue phase, that have true three- 
dimensional periodic order are mcsophases. 

7. nematic: a liquid crystalline material composed 
of anisotropic (rod or disk shaped) particles with 
long-range orientational but no long-range transla- 
tional order. This term comes from the Greek word 
veiMjoa for thread. A nematic is often filled with de- 
fects that look like threads under cross polarizers. 



nematogen: 

phase. 



a molecule that forms a nematic 



9. smectic: from the Greek word afiej/ia for soap. 
A smectic phase is a "solid" in one dimension and 
a fluid in the other two directions. It consists of 
equally spaced parallel layers. 

10. steric: arising from hard-core, excluded- volume in- 
teractions. This term comes from the Greek word 
arepeoa for solid. 



APPENDIX A: EXPRESSION FOR THE TORQUE FIELD 

In this appendix we recast the expression for h in a 
simplifying limit, viz. when the molecules are perfectly 



aligned along the local nematic axis (but their biaxial axis 
is not fixed). We start from Eq. (34). As mentioned in 
the text, it is clear that h is zero if the system is achiral. 
Thus we are interested in the terms in the free energy 
which arc linear in the chiral parameters, ip n . So we may 
write 



h 



d 2 F 



dqdip n 



(Al) 



Note that q enters the calculation in the following way. 
We consider a helical phase in which the director n is 
given by Eq. (1). In the small q limit, the atomic co- 
ordinates are displaced by an amount 5r\ a from their 
reference positions in the nematic phase, where Sr\ a = 



e ljk Suj 



5uj, 



qe J (e • R A ) and where e is an arbi- 



trary unit vector perpendicular to n. In this sense the 
total potential energy U has a g-dependence such that 



dU 

dq 



y dU dr\ a 



(A2) 



With this understanding, one evaluates Eq. (Al) as 



4e^» 



/ d 2 u 

\ dqdip n 



+ Xn 



(A3) 



where (•) denotes a thermodynamic average in which 
the density matrix cxp[— U/{k B T)\ is evaluated when all 
molecular chirality is turned off and 



1 



Xr 



dU dU 



k B T \ dq dip n 



(A4) 



Note the appearance in h of the terms in \ n . As we 
will show elsewhere (Harris, et al. 1999), these terms, 
which normally arc not considered, are needed to obtain 
the expected result that h vanishes in the limit of an 
isotropic fluid (for which the nematic order parameter 
vanishes). However, in the limit of nearly complete ne- 
matic order (which we consider here), these terms in \ n 
are negligible. Superficially it may appear that we have 
to isolate the dependence of U on the chiral parameters. 
However, since achiral components of U do not contribute 
to (dU/dq), this step is, in fact, not necessary, so that 



(A5) 



This equality is the basis of our calculation - it allows us 
to calculate OF/ dq microscopically in terms of molecular 
interactions. In particular, note that the result is ex- 
pressed in terms of a correlation function to be evaluated 
in the nematic (q = 0) limit (which we do implicitly in 
the following.) It is important to note that in the more 
realistic limit when the orientations of the molecules fluc- 
tuate away from the local nematic direction, both terms 
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in Eq. (A3) must be retained. Such a calculation has not 
yet been carried out. 

Writing the potential energy as a sum of identical (for 
simplicity), pairwise central- force interactions V^(R) and 
using Eqs. (33), (A2), and (A5), we find 

qo = ~ 2K T \ ^ € iokdiV {H B +r B[j -'R A -r Aa ) 

2 \ BA0OL 

ej -{(e-R fl )r*0-(e-R A )r* Q }^ . (A6) 

Because the system is uniaxial, we can average over all 
perpendicular directions e, so that e^e^ — ► \ — niUj). 
In this case Eq. (A6) becomes 




where Rj^ is the projection of R onto the plane per- 
pendicular to n and R = R B — R^. By translational 
invariance, the second of the two terms inside the square 
brackets of Eq. (A7) must vanish: for fixed molecule B, 
one can shift the origin of the coordinate system by a 
fixed vector A. This would make the second term de- 
pend on the choice of origin, which it cannot. Hence we 
find that 




where t ba is the torque exerted on molecule B by 
molecule A: 

t 'ba = e i]kr 3 B(j d k V (R + r Bfj - r Aa ) . (A9) 

Pa 
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